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VAPOR  PRESSURE  OF  VX 


1.  INTRODUCTION 

The  vapor  pressure  of  VX  has  been  measured  previously  between  108  and  231°C  using 
differential  thermal  analysis^  and  between  52  and  100°C  using  Knudsen  effusion^.  The  lower 
limit  of  previous  work  was  52°C  where  the  vapor  pressure  was  reported  to  be  0.01 1 8  Torr 
(1.57  Pa).  While  equations  are  available  to  extrapolate  beyond  the  experimental  range,  the  errors 
associated  with  these  calculated  vapor  pressures  increase  as  the  extrapolation  gets  farther  from 
measured  data.  Therefore,  it  is  desirable  to  have  measured  data  over  as  wide  a  range  as  possible. 
The  current  work  employs  a  modified  ASTM  gas  saturation  method^  and  extends  the  temperature 
range  for  measured  VX  vapor  pressure  data  to  -13°C,  which  allows  vapor  pressure  measurements 
for  VX  well  over  3  orders  of  magnitude  lower  than  previously  reported. 

Vapor  pressure  data  are  used  for  a  wide  variety  of  applications  including  assessing 
evaporation  rates  at  ambient  conditions,  calibration  of  detectors,  vapor  generator  design  for 
evaluation  of  physiological  effects  resulting  from  exposure  to  toxic  vapors,  and  evaluation  of 
vapor  filtration  performance.  One  method  used  to  generate  and  control  the  concentration  of 
vapors  of  interest  in  a  gas  stream  is  to  pass  an  inert  carrier  gas,  e.g.,  dry  nitrogen  or  helium,  over 
a  condensed  material  maintained  in  a  container  at  a  controlled  temperature.  If  the  ambient 
pressure,  carrier  flow  rate,  and  vapor  pressure  of  the  liquid  are  known,  the  concentration  exiting 
the  generator  can  be  calculated  by  use  of  the  ideal  gas  law.  Conversely,  if  the  ambient  pressure 
and  carrier  gas  flow  rate  are  known  and  tiie  concentration  can  be  measured,  the  vapor  pressure 
can  be  inferred.  Once  the  vapor  pressure  dependence  on  temperature  is  fully  characterized, 
effluent  concentration  can  be  precisely  controlled  by  changing  the  carrier  flow  rate  and 
temperature.  Accordingly,  it  is  important  to  have  an  accurate  knowledge  of  the  temperature- 
vapor  pressure  relationship  for  each  vapor  of  interest. 


2.  EXPERIMENTATION 

Saturated  VX  vapor  streams  were  generated  by  flowing  nitrogen  carrier  gas  at  30.0 
standard  cubic  centimeters  per  min  (seem)  through  a  glass  vessel,  i.e.,  saturator,  containing  liquid 
VX,  shown  schematically  in  Figure  1.  The  saturator  has  two  vertical  7-mm  o.d.  tubes, 
about  20  cm  long,  connected  to  either  end  of  a  100-mm  long,  25-mm  o.d.  cylindrical  horizontal 
tube,  which  held  the  liquid.  The  horizontal  tube  contained  a  concentric  hollow  ceramic  cylinder 
in  contact  with  the  liquid  sample  fabricated  in  a  manner  to  allow  the  carrier  to  make  three  passes 
along  the  length  of  the  tube  prior  to  exiting  and  to  increase  the  contact  area  since  the  sample 
wetted  the  rough  ceramic  surface.  This  geometry  provided  for  a  carrier  gas  residence  time  of 
about  1  min  in  the  saturator. 

In  this  work,  the  saturator  was  loaded  with  7  g  of  neat  VX,  Chemical  Agent  Standard 
Analytical  Reference  Material  (CASARM),  Lot  Number  VX-U-2128-CTF-N,  94.6%,  used  as 
received.  The  VX  sample  used  for  these  measurements  was  analyzed  before  data  were  collected 
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by  NMR  spectroscopy,  Gas  Chromatography  (GC)-Mass  Spectrometry  (MS),  and  GC  using  a 
Thermal  Conductivity  Detector  (TCD).  The  VX  sample  used  in  the  vapor  pressure  study  was 
also  analyzed  immediately  after  the  vapor  pressure  data  were  collected  by  GC/TCD  with  no 
observable  change  in  purity. 

The  saturated  VX  stream  was  sampled  by  drawing  5-20  seem  of  the  saturator  effluent  to 
the  Dynatherm  concentrator,  a  modified  ACEM  Model  900  (Dynatherm  Analytical  Instruments 
Inc.,  Kelton,  PA)  system  containing  a  tenax  adsorbent  tube.  The  flow  rate  was  controlled  using  a 
calibrated  10-  or  100-sccm  mass  flow  controller  (Tylan  Model  FC-280,  Austin,  TX),  which  was 
connected  to  the  house  vacuum  for  1-25  min,  depending  on  the  amoimt  of  VX  to  be  collected, 
i.e.,  flow  rate  and  temperature.  The  saturator  effluent  was  sampled  at  a  "tee"  connection  3  cm 
downstream  of  the  saturator  constant-temperature  bath.  The  sampling  line  was  constructed  of 
1/16-in.  o.d.,  0.040-in.  i.d.  Silico  Steel  (fused  silica  lined  stainless  steel,  Restek  Corporation., 
Bellefonte,  PA).  All  connectors  between  the  saturator  effluent  line  "tee"  and  the  concentrator 
were  316  stainless  steel.  The  concentrator  sampling  line  between  the  "tee"  and  the  ACEM  was 
heated  to  ca.  100°C  to  reduce  sampling  line  losses  and  thus  reduce  the  time  needed  to  achieve 
steady  state  following  sample  temperature  changes.  The  ACEM  900  valve  box,  external  six-port 
valve,  and  transfer  line  to  the  GC  were  maintained  at  150®C.  The  sample  was  concentrated  by 
the  ACEM  900  using  the  following  operating  parameters.  The  temperature  of  the  10-mm  o.d. 
tenax  collection  tube  was  maintained  at  35°C  during  the  variable  1-25-min  sample  collection 
time.  A  1-min  purge  using  dry  nitrogen  was  employed  as  the  next  step  of  the  sampling  process. 
The  sample  was  rapidly  heated  to  275°C  under  a  flow  rate  of  20  seem  UHP  grade  heliiun  for 
3  min  and  transferred  to  the  ACEM  900  tenax  focusing  trap,  which  was  maintained  at  35°C. 
Transfer  continued  for  an  additional  min  as  the  10-mm  o.d.  tenax  collection  tube  cooled.  The 
focusing  trap  was  then  rapidly  heated  to  300°C  under  a  flow  of  8  seem  UHP  grade  helium  for 
3  min  to  effect  sample  transfer  to  the  GC  column.  The  GC  column  (Restek  Coiporation., 
Bellefonte,  PA),  a  30-m  x  0.53-mm  i.d.  fused  silica  column  with  1.0-pm  RTx-1  (polydimethyl- 
siloxane),  was  maintained  at  50'’C  for  2  min,  then  heated  at  a  rate  of  10°C/min  to  200°C,  and 
maintained  at  200°C  for  5  min  well  after  the  VX  had  eluted.  A  Hewlett  Packard  (HP)  Model 
5890  Series  II  GC  equipped  with  a  flame-photometric  detector  and  a  phosphorus-selective  filter 
(FPDP)  was  used  as  the  analytical  detection  system  in  this  work.  Helium  was  used  as  the  GC 
carrier  gas  at  a  flow  rate  of  8  seem  and  as  detector  make-up  gas  at  a  flow  rate  of  12  seem. 
Combustion  gases  used  to  support  the  FPDP  were  air  at  1 00  seem  and  hydrogen  at  75  seem. 

Using  the  instrumentation  and  operating  conditions  described,  VX  eluted  at  15.5  min,  which 
represents  a  GC  colmrm  temperature  of  185°C.  The  FPDP  minimum  detectable  limit  for  VX  was 
determined  to  be  0.3  hg  with  a  signal-to-noise  ratio  >2.  Initial  saturator  GC  analyses  revealed 
more  than  20  phosphorus-containing  compounds.  Surprisingly,  the  VX  peak  measured  from  the 
saturator  effluent  inittially  accounted  for  only  about  7%  of  the  total  GC  peak  area  as  opposed  to 
the  calibration  and  liquid  analysis  data  which  showed  that  the  VX  peak  accounted  for  95%  of  the 
total  area.  The  prominence  of  the  VX  GC  peak  fi:om  the  saturator  effluent  increased  slightly 
during  the  2  weeks  of  data  measurement  to  a  maximum  of  9%  of  the  total  area.  The  low 
abundance  of  VX  in  the  vapor  phase  is  a  consequence  of  its  low  volatility  compared  to  the 
known  impurities,  most  of  which  have  considerably  higher  volatility  than  VX.  The  slow  clean 
up  of  VX  has  also  been  observed  in  our  laboratory  for  similar  low-volatility  materials.  This 
phenomenon  does  not  present  a  problem  for  vapor  pressure  calculation  as  discussed  below. 
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A  “stock”  solution  for  calibration  was  prepared  by  adding  5  pL  of  neat  VX  to  about 
20  mL  of  isopropyl  alcohol  (IP A)  using  a  5-pL  Drummond  Model  105  Microdispenser 
(Drummond  Scientific  Company,  Broomall,  PA)  and  diluting  to  25  mL  with  IPA.  Converting 
this  volume  of  VX  to  mass  using  the  liquid  density  at  room  temperature  (1.0103  g/ml  at  22.5°C) 
and  correcting  for  VX  purity  (94.6%)  resulted  in  a  solution  with  a  concentration  of  1 91 .2  ng  VX 
per  microliter  of  IPA.  A  calibration  standard  was  prepared  from  this  “stock”  solution  to  span  the 
experimental  data  range  of  interest.  A  7.648-ng/pL  (28.6  micromolar)  standard  was  prepared  by 
diluting  400  pL  of  the  “stock”  solution  to  10  mL  with  IPA  using  a  200-pL  (Drummond  Model 
385)  Microdispenser  (Drummond  Scientific  Company,  Broomall,  PA).  The  “stock”  solution  and 
calibration  standard  were  prepared  using  volumetric  flasks  that  had  been  preconditioned  with 
similar  concentrations  of  VX  in  IPA  for  about  2  weeks. 

All  calibrations  were  performed  using  standards  prepared  from  the  stock  solution  on  the 
same  day.  The  FPDP  was  calibrated  by  making  1-6-pL  injections  of  VX-isopropyl  alcohol 
standards  onto  the  lO-mm  o.d.  tenax  collection  tube  mounted  horizontally  in  the  Dynatherm 
Model  10  Tube  Conditioner  (Dynatherm  Analytical  Instruments,  Incorporated,  Kelton,  PA)  at  a 
laboratory  temperature  of  22.5“C  under  a  30-sccm  flow  of  nitrogen  carrier  gas  for  2  min.  After 
each  injection  was  made,  the  tenax  tube  was  placed  into  the  ACEM  900  for  transfer  to  the  GC 
column  for  analysis.  ACEM  900  and  GC  operating  parameters  were  identical  to  those  used  for 
experimental  data  acquisition.  The  resulting  calibration  curve  was  generated  by  plotting  FPDP 
area  counts  versus  mass  of  VX  injected  and  was  slightly  non-linear.  The  quadratic  equation 
shown  in  Equation  1  described  the  calibration  data  most  accurately. 

Equation  (1)  A  =  0.003045»mvx^  +  0.3875»mvx 

where 

mvx  =  mass  ofVX(ng)  injected 

A  =  GC  area  (millions  of  area  counts) 

The  resulting  raw  data  and  calibration  curve  are  shown  in  Figure  2.  The  slight 
nonlinearity  is  due  to  the  presence  of  sulfur  in  VX.  Similar  calibrations  performed  in  our 
laboratory  for  compounds  containing  both  phosphorus  and  sulfur  have  produced  similar  results, 
while  calibrations  for  compounds  containing  phosphorus  only  were  more  linear.  The  reason  for 
this  behavior  is  overlap  of  the  low-energy  tail  of  the  sulfur  molecular  emissions  with  the 
phosphorus  filter  window  near  530  nm.'*  As  seen  in  Figure  2,  the  data  scatter  of  the  calibration 
curve  was  minimal  (R^  =  0.999).  Two  VX  calibrations  performed  several  months  apart  produced 
similar  calibration  results.  The  linear  and  quadratic  coefficients  of  the  calibration  curve  are 
indicative  of  the  selectivity  of  the  phosphorus  filter.  In  this  case,  the  phosphorus  selectivity 
appears  to  be  higher  than  99%,  although  this  fraction  will  also  depend  on  other  factors,  including 
relative  intensity  of  the  phosphorus  and  sulfur  emissions  and  specific  filter  characteristics. 

The  temperature  of  the  sample  container  was  controlled  by  immersing  it  in  a  water- 
ethylene  glycol  bath.  The  temperature  of  the  bath  was  measured  to  within  0. 1  °C  using  a 
calibrated  thermometer.  The  ambient  pressure  was  measured  periodically  during  each  run  using 
a  Princo  Instruments  Nova  Model  mercury  manometer  (Princo  Instruments,  Incorporated, 
Southampton,  PA).  All  manometer  readings  were  corrected  for  temperature  and  latitude 
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according  to  the  manufacturer's  directions.  These  readings  were  used  in  the  vapor  pressure 
calculation  as  described  below.  Owing  to  the  experimental  design  and  operating  parameters, 
there  should  be  very  little  over-pressure  inside  the  saturator.  This  correction  is  believed  to  be 
very  small  and  has  not  been  used  in  this  work. 

Calculation  of  VX  vapor  pressure  from  the  mass  collected  at  each  temperature  was 
performed  using  the  ideal  gas  law.  The  number  of  moles  of  nitrogen  exiting  the  saturator  was 
determined  from  the  concentrator  flow  rate  (seem)  and  the  sampling  time.  VX  mass  loss 
measured  in  this  study  was  kept  in  the  30-50  ng  range  by  using  longer  sample  collection  times 
and  higher  flow  rates  for  the  lower  temperature  measxirements.  The  VX  mass  was  determined  by 
comparing  GC  areas  to  those  measured  in  the  calibration  studies.  The  derived  mass  was  used  to 
determine  the  number  of  moles  of  VX  evaporated.  The  total  number  of  moles  exiting  the 
saturator,  i.e.,  VX  and  carrier  gas,  is  the  sum  of  the  two.  VX  vapor  pressure  was  calculated  using 
the  ideal  gas  law  assumption  as  shown  in  equation  (2).  The  calculated  vapor  pressure  was 
corrected  for  the  sample  purity  by  taking  Raoult's  Law  into  account. 


Equation  (2)  Pvx  =  (Pamwent  *  (nvx/(nvx  +  n„^„)))/yvx 


where 


Pvx  ~ 

P  = 

^  ambient 

iivx  ~ 

Pearlier 

P 

R 

T 


Yvx 


VX  vapor  pressme  in  pascals 

corrected  barometric  pressure  reading  (Pa) 

moles  VX  =  measured  VX  mass  loss  (g)/(267.4  g/gm  mole) 

moles  carrier  =  calibrated  flow  rate  (seem)  •  time  (minutes)  •  (P/RT) 

101325  Pa 

8314000  cc-Pa/mole-°K 
273.16  °K 

VX  mole  fraction  (0.946) 


3.  RESULTS 


Analyses  of  the  CASARM  sample  performed  prior  to  data  collection  indicated  a  94.6% 
mole  fraction  of  VX.  Most  of  the  significant  impurities  had  been  identified  and  are  common  to 
VX  samples.  Most  of  these  impurities  have  considerably  higher  vapor  pressure  than  VX  and 
were  reduced  by  prolonged  purging  with  the  carrier  gas.  The  sample  was  purged  for  several  days 
at  room  temperature  before  initiation  of  data  collection  to  reduce  the  concentration  of  the 
impmities.  This  effort  was  partially  successful  as  indicated  by  the  relative  decrease  in  the  GC 
area  attributable  to  impurities.  A  much  higher  success  could  have  been  achieved  for  a  smaller 
sample  size,  and  future  saturator  cell  design  will  incorporate  changes  to  allow  smaller  samples  to 
be  used. 

Vapor  pressure  measurements  were  performed  using  the  automatic  cycling  capability  of 
the  vapor  concentrator  and  sampling  system.  Typically,  each  sample  was  collected  for  between 
1  and  25  min  at  5-20  seem  onto  the  tenax  adsorbent.  Following  sample  collection  and  transfer  to 
the  focusing  trap,  the  adsorbate  was  transferred  to  the  GC,  and  the  GC  column  was  heated  from 
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50°C  to  200°C  over  a  22-min  period.  This  process  took  almost  1  hr,  and  the  process  was 
repeated  vmtil  the  area  under  the  GC  peak  for  VX  was  stable.  Usually,  it  took  5-8  hr  for  the  VX 
area  to  achieve  stability  following  sample  temperature  changes,  while  the  higher  volatility 
components  responded  and  stabilized  much  more  quickly. 

The  raw  and  corrected  vapor  pressure  data  measured  in  this  work  are  shovm  in  Table  1 . 
Each  entry  represents  an  average  of  at  least  10  separate  measurements.  Since  the  liquid  VX 
sample  mole  fraction  was  only  about  95%,  the  measured  vapor  pressures  were  corrected  using 
Raoult's  Law.  The  majority  of  the  data  reported  here  were  measured  at  a  saturator  flow  rate  of 
30  seem.  Two  experiments  were  performed  at  50  seem  saturator  flow  rate  once  the  GC  area 
counts  had  stabilized  for  the  20°C  data  point.  The  insensitivity  of  calculated  VX  vapor  pressure 
to  flow  rate  indicates  achievement  of  saturation,  even  at  the  higher  flow  rate. 

The  VX  vapor  pressure  data  measured  using  the  saturator  in  this  work,  along  with  data 
presented  by  Savage  and  Fielder  and  by  Belkin  and  Brown  are  plotted  in  Figure  3.  A  new 
Antoine  correlation  was  derived  based  on  all  data  considered  in  this  report.  Comparison  of  the 
measured  data  and  calculated  vapor  pressures  over  the  range  imder  consideration  are  shown  in 
Figure  3.  New  Antoine  coefficients  are  listed  in  Table  2  along  with  comparisons  between  the 
measmed  and  calculated  data  based  on  the  new  Antoine  analysis.  Table  3  presents  calculated 
VX  vapor  pressure  based  on  the  previous^  and  new  Antoine  analyses  and  calculated  VX  volatility 
over  the  ambient  temperature  range. 

4.  DISCUSSION 

Previous  vapor  pressure  measurements  for  dimethyl  methylphosphonate  (DMMP)  and 
distilled  water  were  performed  in  our  laboratory  in  a  similar  fashion  by  directly  measuring  mass 
loss.^  This  simple  method  worked  well  for  those  compounds  since  their  vapor  pressures  are  in  a 
range  where  sufficient  weight  loss  was  observed  over  a  reasonable  period  of  time.  The  water 
vapor  pressure  results  were  in  excellent  agreement  with  literature®  data.  The  DMMP  data  were 
consistent  with  earlier  data  and  lowered  the  measured  data  range  by  about  80°C.’  The  direct 
mass  loss  method  is  not  suitable  for  compounds  with  vapor  pressure  in  the  range  of  VX  in  the 
temperature  range  of  interest  since  measurable  changes  would  require  either  weeks  or  months 
instead  of  hours  or  days.  Another  serious  problem  for  weight-loss  methodology  for  compounds 
such  as  VX  is  the  presence  of  higher-volatility  impurities.  As  a  result,  the  current  measurements 
have  resorted  to  measuring  mass  loss  indirectly  using  a  much  more  sensitive  purge-and-trap  gas 
method  coupled  with  gas  chromatographic  analysis  v^ith  FPDP  detection  to  separate  VX  from  its 
impurities. 

Prior  vapor  pressure  measurements  performed  in  this  laboratory  for  DMMP  used  two 
saturator  vessels  very  similar  to  the  one  used  in  this  work.  Those  were  placed  in  series  to 
lengthen  the  contact  time  between  carrier  and  liquid  in  an  attempt  to  ensure  vapor-liquid 
equilibrium.  That  work  Was  performed  using  higher  flow  rates  than  those  used  in  the  present 
work.  Measurements  in  the  earlier  work  showed  that  vapor-liquid  equilibrium  was  achieved  and 
that  very  little  DMMP  was  lost  from  the  second  saturator  for  flows  as  high  as  80  seem,  which 
represented  the  upper  flow  rate  limit  of  the  earlier  work.  For  a  variety  of  reasons,  including  the 
more  sensitive  analytical  method  and  the  hazards  associated  with  using  large  quantities  of  VX,  a 
single  saturator  was  used  in  the  present  work.  To  demonstrate  that  the  vapor  exiting  the  saturator 
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was  in  equilibrium  with  the  liquid,  several  VX  vapor  pressure  measurements  were  done  at 
saturator  flow  rates  of  30  and  50  seem.  For  these  measurements,  the  derived  vapor  pressure  was 
independent  of  flow  rate,  demonstrating  that  equilibrium  had  been  aehieved  at  flow  rates  in  this 
range  using  a  single  saturator. 

The  saturation  method  used  to  measure  VX  vapor  pressure  in  this  work  has  a  very 
significant  advantage  over  mass-loss  methods  such  as  direct  weighing  and  Knudsen  effusion. 
That  advantage  accrues  from  the  GC  analysis  performed  on  the  effluent.  The  GC  analysis  allows 
a  direct  assessment  of  the  analyte  concentration.  In  the  present  case,  VX  is  a  minor  component 
of  the  effluent  stream  owing  to  its  low  vapor  pressure  compared  to  the  myriad  of  impurities  it 
contains.  A  mass-loss  method  can  be  used  for  low-volatility  vapors  with  high-volatility 
impurities  such  as  for  VX,  but  only  with  caution  to  ensure  that  the  sample  is  purified  to  the  point 
that  the  vast  majority  of  the  analyzed  material  is  actually  the  vapor  of  interest.  Analyzing  the 
effluent  and  separating  the  analyte  of  interest  gets  around  this  problem  if  its  mole  fraction  is 
known  since  the  measured  vapor  pressure  can  be  corrected  to  that  of  the  pure  compound  by  use 
of  Raoult's  Law.  The  Antoine  fit  of  the  corrected  data  has  been  derived  for  all  data  listed  in 
Table  1  and  is  presented  as  Equation  3.* 

Equation  (3)  In  VP  (Pascal)  =  a-b/(c+T) 

where  a,  b,  and  c  are  fit  parameters,  and  T  is  the  absolute  temperature. 

The  reduction  of  volatile  impurities  that  accompanied  prolonged  sample  purging  in  the 
saturator  raised  questions  regarding  the  sensitivity  to  this  effect  of  the  methods  used  previously 
to  measure  VX  vapor  pressure.  Differential  thermal  analysis  (DTA)  measurements  are  carried 
out  under  sufficiently  different  conditions  that  purging  of  volatile  compounds  should  not  affect 
the  results.  However,  the  conditions  for  Knudsen  effusion  measurements  are  similar  enough  to 
the  saturator  experiments  that  this  issue  needs  to  be  considered. 

In  DTA,  the  sample  and  a  thermocouple  are  loaded  into  a  glass  capillary  tube,  which  is 
then  placed  in  the  instrument  heating  block.  The  system  is  stabilized  at  the  desired  pressure,  and 
the  block  is  heated  until  boiling  of  the  liquid  is  detected  based  on  the  time-temperature  plot.  The 
technique  is  rapid,  and  a  data  point  can  be  generated  in  a  matter  of  min.  Pure  materials  produce  a 
sharp  endothermic  peak  with  a  vertictil  leading  edge  as  boiling  occurs  at  a  constant  temperature. 
The  boiling  point  is  defined  as  the  peak  onset  temperature  and  is  taken  at  the  intersection  of 
tangents  to  this  curve.  The  presence  of  impurities  tends  to  broaden  the  boiling  peaks  as  the 
boiling  range  increases,  thus  changing  the  peak  onset  temperature.  Due  to  the  speed  of  the 
measurement  and  the  fact  that  the  result  depends  on  the  boiling  temperature  of  the  liquid, 
concentration  of  volatile  impurities  in  the  vapor  should  not  impact  vapor  pressure  measurements 
by  DTA.  As  seen  in  Figure  3,  the  low  end  of  the  DTA  data  deviate  significantly  from  the 
Antoine  fit.  It  is  possible  that  slow  pumping  rates  are  responsible  for  the  positive  deviation.  In 
that  case, -the  local  pressure  would  be  slightly  higher  than  the  measured  pressure  and  lead  to 
apparently  higher  than  actual  vapor  pressures.  If  the  DTA  data  below  1000  Pa  are  eliminated, 
the  average  error  for  the  remaining  data  is  reduced  to  below  6%. 

In  the  Knudsen  effusion  method,  vapor  pressure  is  calculated  based  on  weight  loss  over 
time  as  the  test  material  effuses  through  an  orifice  of  known  area  at  a  controlled  temperature 
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under  high  vacuum.  Measurements  are  time-consuming  since  it  can  take  hours  to  achieve 
thermal  equilibrium  under  high  vacuum.  The  material  effusing  from  the  sample  cup  is  not 
analyzed,  so  it  is  difficult  to  determine  how  much  of  the  weight  loss  is  due  to  impurities  rather 
than  the  material  of  interest.  On  the  other  hand,  the  samples  are  small,  so  sample  purification  by 
distillation  should  be  rapid.  The  presence  of  impurities  in  a  sample  is  indicated  by  a  change  in 
the  slope  of  the  weight-vs-time  curve.  Weight  loss  data  are  not  used  for  calculation  of  the  actual 
vapor  pressure  until  a  stabilized  rate  of  weight  loss  has  been  established.  At  this  point,  it  is 
assumed  that  the  effusing  species  is  the  compoimd  of  interest.  However,  this  may  not  be  the  case 
if  the  impurities  are  close  in  volatility  to  the  test  material  or  if  the  sample  contains  a  significant 
quantity  of  a  lower-volatility  component.  Calculation  of  vapor  pressure  under  these  conditions 
assumes  that  the  molecular  weight  of  the  effusing  species  is  known.  If  weight  loss  attributed 
solely  to  the  test  material  actually  includes  impurities,  the  calculated  result  will  be  inaccurate. 

The  Antoine  equation  previously  considered  to  reflect  the  most  accurate  vapor  pressure 
data  for  VX  was  presented  by  Savage  and  Fielder.^  That  report  cites  three  discrete  data  sets  from 
which  that  Antoine  equation  was  generated: 

•  DTA  data,  temperature  range  1 08  to  23 1  “C 

•  Knudsen  data,  temperature  range  52  to  100°C 

•  Knudsen  data,  temperature  range  30  to  40°C 

The  DTA  data  set  attributed  to  Belkin  and  Brown’  will  not  be  discussed  further  here. 
Savage  and  Fielder^  cited  five  laboratory  notebooks  containing  the  experimental  data.  For  the 
higher  temperature  Knudsen  data  set,  only  a  typed  table  of  the  same  temperature/pressure  pairs 
listed  in  the  report  could  be  located  in  the  cited  notebooks.  None  of  the  detailed  raw  data  typical 
of  a  Knudsen  notebook  entry  was  found  for  VX.  Savage’s  notebooks  did  contain  numerous  other 
lengthy  entries  describing  Knudsen  experiments  with  other  compounds.  These  contained 
comments  indicating  awareness  of  the  implications  of  a  changing  slope  with  Knudsen  weight 
loss  data.  This  suggests  that  the  VX  Knudsen  data  were  treated  in  a  similar  manner. 

The  sowce  of  the  second  data  set  was  not  cited  by  Savage  and  Fielder.  Data  between 
30  and  40°C  were  measured  by  Haskell  in  the  mid  1960s  and  were  documented.’  Haskell’s 
notebook  contains  a  discussion  by  Coulter  regarding  which  experimental  points  should  be  used 
for  calculation  of  the  Antoine  equation  and  the  reasons  supporting  his  selections.  This  discussion 
clearly  indicates  that  for  the  VX  measurements  in  this  lower  range,  the  initial  weight  loss  was 
due  to  loss  of  volatile  impurities,  and  those  points  were  discarded. 

The  unavailability  of  raw  data  for  Savage’s  and  Fielder’s  Knudsen  measurements 
complicates  assessment  of  the  literature  data  but  does  not  justify  excluding  those  points  from  a 
new  Antoine  fit  for  the  combined  literature  and  saturator  data  sets.  If  the  temperature  range  of 
Savage’s  and  Fielder’s  data  covered  the  high  or  low  end  of  the  temperature  range  being 
considered  here,  an  argument  could  be  made  to  exclude  those  points.  Since  the  points  fall  in  the 
middle  of  the  experimental  range  of  all  data  and  have  good  overall  consistency  with  the  other 
data  sets,  they  were  included  in  calculation  of  a  new  VX  Antoine  equation. 

The  new  Antoine  fit  to  the  saturator  data  produced  an  average  difference  between 
calculated  and  measured  vapor  pressures  of  1 .8%  for  the  saturator  data,  7.0%  for  the  Knudsen 
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data,  and  1 1 .6%  for  the  DTA  data,  with  an  overall  average  difference  of  9.0%.  The  Antoine 
equation  developed  by  Savage  and  Fielder  is  lower  by  a  factor  of  three  than  the  lowest 
temperature  saturator  data  to  over  1  order  of  magnitude  lower  at  -40°C  (estimated)  as  seen  in 
Tables. 

5.  CONCLUSIONS 

This  work  represents  the  first  time  that  the  vapor  saturation  method  in  combination  with 
purge-and-trap  chromatography  has  been  used  to  measure  the  vapor  pressure  of  a  CW  agent.  The 
methodology  employed  has  a  number  of  advantages  over  methods  previously  used  for  measuring 
the  vapor  pressure  of  toxic,  low- volatility  materials.  The  vapor  saturation  method  in  combination 
with  very  sensitive  analytical  methods  such  as  the  ones  used  in  this  study  shows  great  promise 
for  vapor  pressure  measurements  well  below  10'^  Pascal.  In  addition,  the  analytical  methodology 
used  in  the  present  work  gets  around  sample  and  effluent  purity  problems  since  the  vapor(s)  of 
interest  are  separated  by  GC  and  analyzed  directly.  The  vapor  pressure  of  the  pure  compound 
can  be  derived  using  Raoult's  Law  if  the  liquid  sample  purity  is  known. 

The  saturator  vapor  pressure  data  measurements  have  lowered  the  range  of  measured  VX 
data  by  65°C.  Data  measured  in  this  work  between  -13  and  20“C  indicate  a  higher  (ca.  33%)  VX 
volatility  at  25°C  than  predicted  by  the  most  recent  previous  analysis.^ 

Results  of  this  type  are  important  for  several  applications,  including: 

•  Accurate  assessment  of  liquid  VX  evaporation  rates  at  ambient  conditions 

•  Design  of  vapor  generators  for  low-concentration  toxicity  studies 

•  Calibration  of  sensitive  CW  agent  vapor  detectors 

•  Evaluation  of  CW  agent  vapor  filter  efficiency 
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Figure  2.  VX  Calibration  Curve  and  Analytical  Equation 
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Figure  3.  VX  Vapor  Pressure  Plot  Showing  Previous  Literature  Data,  New  Saturator  Data,  and 
Antoine  Fit  to  all  Data 


Table  1.  Measured  and  Corrected  Vapor  Pressure  for  VX  and  Sample  Volumes 


Temperature 

(°C) 

Measured 

VXVP 

(milliPascal) 

Corrected 

VXVP* 

(milliPascal) 

Sample 
Volume  (see) 

20 

62.2 

65.8 

5 

15 

35.1 

37.0 

10 

5 

10.5 

11.2 

40 

-5 

2.75 

2.9 

120 

-13 

0.81 

0.86 

500 

*  Corrected  using  Raoult's  Law  by  dividing  measmed  vapor  pressure  by  VX  mole  fraction 
(0.946). 
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Table  2.  Measured  VX  Vapor  Pressure,  Calculated  Vapor  Pressure  Based  on  Updated 
Antoine  Parameters,  and  Percent  Difference. 


Antoine  Coefficients 


a 

23.72 

b 

6154.9 

c 

-60.16 

Average  Difference 

9.0% 

Calculated  Boiling  Point 

292°C 

TCO 

Measured  VP  (Pa) 

Calculated  VP  (Pa) 

Difference  (%) 

Saturator  Data 

-13.0 

0.00086 

0.000860 

0.0 

-5.0 

0.0029 

0.00281 

3.1 

5.0 

0.011 

0.0109 

2.1 

15.0 

0.037 

0.0377 

1.4 

20.0 

0.066 

0.0672 

2.2 

Knudsen  Data 

52.0 

1.57 

1.633 

3.8 

56.0 

2.29 

2.306 

0.6 

57.3 

2.40 

2.574 

7.3 

59.6 

3.12 

3.120 

0.0 

62.4 

3.68 

3.925 

6.7 

70.2 

6.73 

7.264 

7.9 

74.5 

8.79 

10.05 

14.4 

80.0 

13.3 

15.03 

13.3 

85.7 

25.8 

22.44 

13.1 

88.5 

27.2 

27.17 

0.2 

90.5 

34.7 

31.08 

10.5 

100.2 

62.4 

58.25 

6.7 

DTA  Data 

108.5 

147 

96.74 

34.0 

109.5 

133 

102.7 

23.0 

125.0 

333 

246.3 

26.1 

126.0 

320 

259.9 

18.8 

128.5 

413 

296.9 

28.2 

132.5 

433 

365.7 

15.6 

138.0 

627 

483.5 

22.8 

147.0 

867 

749.6 

13.5 

148.0 

853 

785.9 

7.9 

160.5 

1626 

1391 

14.5 

■  162.0 

1626 

1485 

8.7 

'  176.5 

2853 

2736 

4.1 

181.0 

2920 

3278 

12.3 

186.0 

3986 

3986 

0.0 

187.0 

3986 

4143 

3.9 

188.5 

3986 

4388 

10.1 

198.0 

6173 

6250 

1.3 

212.0 

9692 

10240 

5.7 

■  219.0  ' 

14000 

12950 

7.5 

221.0 

14000 

13830 

1.2 

221.5 

14000 

14060 

0.4 

231.0 

18120 

19040 

5.1 

231.5 

18120 

19340 

6.7 

237.0 

24490 

22900 

6.5 
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Tables.  Calculated  Vapor  Pressure  of  VX  Over  the  Ambient  Temperature  Range 

Calculated  Using  Antoine  Equations  Derived  from  all  Data  and  by  Savage  and 
Fielder.  Volatility  calculated  based  on  this  v^ork.  (Asterisks  indicate 
extrapolation). 


Temperature 

(°C) 

Vapor  Pressure  (Pa) 
This  Work 

Vapor  Pressure  (Pa) 
Savage  and  Fielder 

Calculated 
Volatility  (mg/m^) 
Based  on  This  Work 

-40* 

7.06E-06 

5.92E-07 

9.73E-04 

-35* 

1.92E-05 

2.19E-06 

2.59E-03 

-30* 

4.93E-05 

7.39E-06 

6.52E-03 

-25* 

1.21E-04 

2.30E-05 

0.0156 

-20* 

2.82E-04 

6.64E-05 

0.0358 

-15* 

6.30E-04 

1.79E-04 

0.0785 

-10 

1.36E-03 

4.55E-04 

0.166 

-5 

2.81E-03 

1.09E-03 

0.337 

0 

5.63E-03 

2.49E-03 

0.662 

5 

0.0109 

5.43E-03 

1.26 

10 

0.0206 

0.0113 

2.34 

15 

0.0377 

0.0228 

4.20 

20 

0.0672 

0.0441 

7.38 

25 

0.117 

0.0825 

12.6 

30 

0.199 

0.150 

21.2 

35 

0.332 

0.264 

34.7 

40 

0.542 

0.454 

55.7 

45 

0.869 

0.760 

87.9 

50 

1.37 

1.24 

136 

55 

2.12 

1.99 

208 

60 

3.22 

3.13 

311 
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